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' Abstract 
00 ' 

QQ ' The computation of the radiation flux related to the Hawking tem- 

. perature of a Schwarzschild Black Hole or another geometric back- 

I ground is still well-known to be fraught with a number of delicate 

' problems. In spherical reduction, as shown by one of the present au- 

. thors (W. K.) with D.V. Vassilevich, the correct black body radiation 

Q I follows when two "basic components" (conformal anomaly and a "dila- 

ton" anomaly) are used as input in the integrated energy-momentum 
■ conservation equation. The main new element in the present work is 

the use of a quite different method, the covariant perturbation theory 
of Barvinsky and Vilkovisky, to establish directly the full effective ac- 
^ ■ tion which determines these basic components. In the derivation of W. 

K. and D.V. Vassilevich the computation of the dilaton anomaly im- 
plied one potentially doubtful intermediate step which can be avoided 
here. Moreover, the present approach also is sensitive to IR (renormal- 
isation) effects. We realize that the effective action naturally leads to 
expectation values in the Boulware vacuum which, making use of the 
conservation equation, suffice for the computation of the Hawking flux 
in other quantum states, in particular for the relevant Unruh state. 
Thus, a rather comprehensive discussion of the effects of (UV and IR) 
renormalisation upon radiation flux and energy density is possible. 
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1 Introduction 



Almost three decades after the (theoretical) discovery of quantum radiation 
from the event horizon of a geometrical background, in particular from a 
Black Hole (BH) ^1 somewhat surprisingly, still the existence of open 
problems is an acknowledged phenomenon - even when only large BH-s are 
considered which, to a good approximation, represent a time-independent 
curved background. 

Actually, the computation of the Hawking effect does not require an anal- 
ysis of the complete evolution of the (massless) fields between infinitely early 
and infinitely late times. It is sufficient to exploit the Energy-Momentum 
(EM) tensor near the future horizon only. The relation between Hawking 
temperature and the radiation flux at infinity is still the object of some de- 
bate. The activity in this field has been rekindled by the work of Bousso 
and Hawking jHi who, on the basis of a computation in 2d dilaton gravity 
resulting from spherically reducing Einstein gravity [3], claimed that an in- 
coming asymptotic flux and thus "anti-evaporation" occurs. Actually already 
some time before ref. [Hj it had been argued [Hj that the so-called conformal 
anomaly (T)2, the trace of the EM tensor in two dimensions^, which had been 
the only input considered^ in ref. 0], cannot provide the complete answer. 

The problem has been focused in refs. where the relation be- 

tween the missing second piece, the 2d "dilaton anomaly" (T^e)^, to the 
pressure component (T^e)^ of the EM tensor in the original 4d theory (in co- 
ordinates time, radius and angles 6, (p) has been established. In the following 
these two essential ingredients of the flux calculation, {T)^ and (T^e)^ (cf. 
eqs. (Pl. lfTIH) below), will be called "basic components" of the EM tensor. 

It had been known for a long time that 2d minimally coupled massless 
scalars - i.e. in the absence of a dilaton field - provide the correct flux from 
the 2d conformal anomaly {T)^ alone. But an actual computation of the 
missing piece (T^q)^ was not available until the work of [T!^ . 

The determination of (T)2 precisely fits into the elegant formalism of heat 
kernel regularisation [T^ I15j. because in 2d "by chance" {T)^ is simply the 
trace anomaly of a massless scalar field which is related to the variation of 
the effective action for a multiplicative (conformal!) factor of the Laplace 
operator. By contrast, the quantum correction (T^q)^ cannot be computed 

lower index 2, 4 attached to an expectation value means computation in spherically 
reduced gravity or directly in four dimensions, respectively. 

^The correct expression for the conformal anomaly in the presence of a 2d dilaton field 
for spherically reduced gravity and for general 2d dilaton theories JH! '7^ taken alone 
yields the same unphysical flux. For a comprehensive review of general 2d dilaton models 
we refer to |S|. 



1 



as easily in this formalism, because it lacks that key property. For this 
reason in refs. the Laplace operator, including a general coupling 

to the dilaton field, has been split into a product of two Dirac operators 
whose combined determinant need be evaluated in fiat space only. Then a 
(multiplicative!) variation with respect to the dilaton field is observed and 
allows the computation of the component (T^g)^, dubbed "dilaton anomaly" 
in that work. Using both basic components in a ( "dilaton-deformed" ) EM 
conservation equation jHl E] produced exactly the black body Hawking fiux 
at infinity from 2d minimally coupled scalars which would follow from the 
Hawking temperature computed from, say, the surface gravity at the horizon. 
The outgoing fiux at infinity precisely coincided with the Stefan Boltzmann 
law for d = 2, which is related easily to the equally correct 4d result. The 
only further input had been the vanishing of the asymptotic incoming fiux 
together with the condition of finiteness (or at least integrability) of the fiux 
at the horizon (Unruh vacuum) in global Kruskal-Szekeres coordinates. In 
addition, by direct functional integration an expression for the total effective 
action in the presence of a general dilaton coupling was obtained. Despite of 
these satisfactory final results, the step sphtting the Laplace operator into 
two linear Dirac operators from a rigorous mathematical point of view seemed 
to be a doubtful one. 

Therefore, in our present paper we attempt to close that loop hole in 
an alternative derivation by a new application of the covariant perturbation 
theory, introduced by Barvinsky and Vilkovisky (BV) pTlllBj . This technique 
allows to proceed directly to the effective action. In ref. [11] it already 
has been shown that the BV effective action reproduces the correct trace 
anomaly of conformally coupled scalar fields in d = 4. Actually, it contains 
more information than just the trace anomaly, although the latter in general 
is the only expectation value that can be calculated directly. In particular, 
the dilaton effective action derived in the present work not only produces the 
correct 2d trace anomaly but also represents another derivation of the dilaton 
anomaly. However, beside the expression found in [TH| ITH] we encounter 
an important IR-renormalisation effect, i.e. something which had been by- 
passed altogether in the only UV-sensitive previous approach. 

That action is used by us only to determine the basic components {T)^ 
and (T^q")^, while the remaining components are reconstructed by integra- 
tion of the EM conservation equation, following the procedure introduced by 
Christensen and Fulling (CF) [T^ . 

In Section 2 we present the dilaton model and some characteristic features 
of spherically reduced (SR) gravity. 

Section 3.1 is devoted to the computation of the effective action of the 
dilaton model by the covariant perturbation theory of refs. |17| I18j. 
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In Section 3.2 we discuss the ambiguity of the non-local effective action 
by a Green function perturbation theory and fix it by appropriate boundary 
conditions and some infrared regularisation. 

The expectation values of the basic components to be derived from the 
effective action are the subject of Section 4 and the remaining components 
in the Unruh state are determined by fixing the constants Q and K of ref. 
[12] accordingly. 

In the Conclusions (Section 5) we summarize and discuss the obtained 
results. 

This paper further contains four Appendices: in Appendix A we shortly 
present the SR procedure and derive some useful formulas. Appendix B con- 
tains the second and third order of the Green function perturbation theory. 
In Appendix C we show the (no n-) conservation of the 2d dilaton EM tensor 
at the quantum level. In Appendix D the regularisation of the heat kernel is 
demonstrated. 

Important basic calculations and concepts of this paper can be found in 
more detail in the PhD thesis ^0] of one of the authors (D. H.), however a 
more careful discussion of (UV and IR) renormalisation issues is presented 
here. 



2 Dilaton Model 



A massless scalar field S is considered on a four-dimensional Schwarzschild 
spacetime M with coordinates — '"^"^ 
not conformally! ) to gravity: 

^2 



M 



X", 6*, v?), in 4d coupled minimally (but 
{VSf 



X 



IGnG 

In the following we set c = G = h = kB = 1. This model can be spherically 
reduced to a dilaton model on a two-dimensional spacetime L by integrating 
out the isometry coordinates 9, ip, (cf. Appendix A 
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On a four- dimensional spherically symmetric spacetime the EM tensor ex- 
hibits only four independent components ^2] = T^e = T'^^) in a vier- 
bein basis 
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where the first block (up to a factor AirX) equals the two-dimensional EM 
tensor T"^ on L. Diffeomorphism invariance in the dilaton model (j^J on-shell 
for S implies the 2d "non-conservation equation" 



dil 



^ 5X 



(4) 



whose solution for the Schwarzschild metric gu = —{grr) ^ = (l ~ ^) 

K 



M2 
1 



(5) 
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(1_2M) [ M2 



2M 
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2M\ drX 5LZ 



dil 
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(6) 



only depends on the integration constants Q, K and the unknown functions 
T, -jY^ to be identified with the basic components. Eq. (gj) is just another 



expression for the 4d conservation equation if one identifies 



1 



6L 



dil 



AttJ^ 6X 



(7) 



This relation can be checked easily for the action Q using the definition 
of the EM tensor T^^ = -^^-^j^- In Appendix C we show that in a fixed 
classical background the non-conservation equation also holds at the level of 
(renormalized) expectation values, 



[m' 



VgX 5W 



(8) 



where W is the generating functional of connected Green functions for the 
dilaton theory which at the one-loop level coincides with the effective action 
(when the propagators of the external lines are amputated). At the quantum 
level the basic components are calculated by variation of the effective action^: 



{T)=9 



5W 



1 5W 
7=f 5^ 



(9) 
(10) 



■^It should be emphasized that is a two-dimensional action and expectation values 
derived from it could differ from those calculated in Ad. 



4 



In the following all quantities are two-dimensional if no dimension index is 
attached. T^t is the flux component of the EM tensor and differs from the 4d 
flux (like all components of the EM tensor) by a factor (47rX)~^ (the dilaton 
X is hidden in the spacetime measure \f--gi = ^\/~9'i)- The constants 
Q, K remain to be fixed by the boundary conditions of the fields and are 
thus related to the quantum state of the system |T2]. 

All physical states are characterized by the choice (5 = which is a 
necessary condition for the finiteness of the EM tensor at the horizon in global 
coordinates. The Hartle-Hawking state is given by vanishing total flux 
i^H = (thermal equilibrium), whereas the Unruh state \U) is determined 
by vanishing incoming flux (leading to a non-zero Ku). The (unphysical) 
Boulware state \B) is defined by vanishing fields in the asymptotic region. 
This is accomplished by setting Kb = and fixing Qb ^ ^ appropriately. 
Although not to be interpreted as a physical state, because of its natural 
boundary conditions, in a certain sense it nonetheless is the natural state of 
the effective action. On the one hand, in order to represent a well-defined 
integral over the fields, natural (vanishing) asymptotic values for them are 
necessary. On the other hand, in the path integral (cf. ()26|) below) the field S' 
is a sum of a classical solution 5*0 and a quantum correction Sg. The standard 
procedure is to set 5*0 = - otherwise one would have surface terms that 
would make the application of the heat kernel method very difficult. This 
means that the incoming and outgoing states correspond to the vacuum, i.e. 
the Boulware state. If the expectation value of the EM tensor is calculated 
from the effective action with the Boulware state values of Kb = and Qb, 
any other quantum state with K = Kb+K, Q = Qb+Q can be reconstructed 
by simply adding (to the first block in 0) a term 



which is a special solution of the (non-)conservation equation. Of course, this 
procedure works only if the basic components are insensitive with respect to 
the state of the effective action which is true, as long as the radiation does 
not affect the spacetime geometry significantly, i.e. in the quasi-static phase 




/ K-Q 
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(11) 



of a BH [20]. 
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3 Non-Local Effective Action 



3.1 Covariant Perturbation Theory 



The relation [TH [1^] between the Eukhdean effective action and the heat 
kernel e"*^"^ for the differential operator O is given by^ 



2 ds 



1 d 1 dr 

-tre 



121 



s=0 



2dsT{s)J, ri- 
The trace of the heat kernel may be expressed in a coordinate basis 

tre-^"= / {x\e-^^\x) ^d^x = [ Go{x,x-T)^d^x. (13) 

J M JM 

In the most common applications the heat kernel is expanded around r = 
[THirni- In contrast to that, the aim of the method developed in pT j IT ^ Hn] 
is to use (fT^ directly in order to find an expression of the heat kernel which 
is valid for all values of r. This allows performing the r-integration and 
computing the effective action for any Euklidean Laplacian 

= -A-E, (14) 

where A = g^'^V^'Vu is the contraction of two general covariant derivatives 
(which may include a gauge part) by some Euklidean metric g, and E is 
an endomorphism, i.e. some linear bounded map from the space of fields 
into itself. The covariant perturbation series is based on a separation of the 
spacetime metric into a flat part g and a perturbing part h: g^^ = g^^ + h^y 
(we use the notation of BV). Nevertheless, each order can be represented by 
covariant expressions corresponding to the full metric g, such as the scalar 
curvature R, because the flat metric does not produce gravitational effects 
R{g) = 0. Then one expands the Laplacian and the heat kernel in orders of 
h: 

O = -Ao + /i^'^V^V, + ..., (15) 

oo 

Go{r) = e-^- = ^GS(r), (16) 

n=0 

where Gq{t) = e^^" and Aq is the fiat Laplacian. Up to the second order in 
the curvature the trace of the heat kernel in even dimensions d = 2uj,uj eN 



We denote Euklidean objects by an index £. 
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is found to be^ \TE\ 
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-Or 



1 



+ T 



R 



(47rr)'- JM 
1 



tr <^ II 



+E 



16(-rA) ■ 32 
/(-rA) , /(-rA) - 1 



I fi-rA) ^ /(-rA)-l ^ 3[/(-rA) - 1] ' ^ 



'-tA] 



8(tA) 



6 



2(-tA) 



f(-rA) f(-rA) 
12 2 



(17) 



where 



In the present paper we restrict ourselves to the second order of covariant 
perturbation theory, i.e. terms up to R^, ER, E"^. This is at least sufficient 
to compute the exact trace anomaly which is completely determined by a 
single term of the local Seeley-DeWitt expansion ^3 ^] corresponding to 
that order. With respect to other expectation values like the dilaton anomaly 
the necessity for higher orders cannot be excluded, a priori, though. 

For d = 2 the trace of the heat kernel (fTTj) (to this order) in (fT^ produces 
five types of integrals. They contain IR or (and) UV divergences that have 
to be regularized by restricting the range of the r-integration as dr in the 
limit T — i> oo,e — >■ 0+. Examples for these rather tedious calculations are 
shown in Appendix D. They can be done analytically to leading order in T, e 
for these IR, resp. UV regularisation parameters, and for the corresponding 
next finite terms. The formal contribution of these terms to the effective 
action (fT^ is given by 
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ds 1 T(s] 
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ds T(s) 



ds { r(S; ^0 
d ( 1 



ds 1 T(s 
1 



r'-^dr 
T'-^dr 



ds 1 r(s; ./Q 



r'f{-TA)dT 
i/(-rA)-l 



-A 



dr 



A2 



'In the following we set lo — 1, i.e. d = 2. 



s=Q 



s=0 



s=0 



s=0 



s=0 



In^ I R,UV (19) 
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e 

2-[ln(-rA)+7g] 

-A 

2 - In(-rA) - 7g 

-A 

ln(-TA)+7^-| 
6(-A) 
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where ~ 0, 57721 is the Euler constant. Collecting all terms in the heat 
kernel p7|) it turns out that the contributions from the divergent parts 
£:~^,lne, and InT surprisingly cancel to zero in those nonlocal terms con- 
taining the scalar curvature R. Thus, the most general two-dimensional 
regularized effective action to second order of perturbation theory in i?, resp. 
E reads 



[9] 



967r 



12co - ci(2i? + 12E) + {R + 12E)^R + 12E^E 



(24) 

Here we have introduced the following regularisation terms: Cq = e ^,Ci = 
ln(£:/T),C2 = ln(— TA) + 7^;. Actually the term oc ciR can be dropped, 
being a total divergence. Note that C2 beside a logarithmic divergence con- 
tains an ill-defined expression in the Laplacian. In Section 4 the (eventual) 
contributions of the regular and divergent terms will be discussed separately. 



3.2 Effective Dilaton Action 

So far all steps in this section were valid for general Laplacian A and endo- 
morphism E. In order to establish the effective action of the dilaton model 
(j21) we have to specify A and E accordingly, while returning to Lorentzian 
spacetime. To achieve this we must reconsider the four-dimensional generat- 
ing functional, determined by the matter part of 



-g4d X 



(25) 

A/" is a normalization constant and the factor -^—g in the path integral mea- 
sure establishes diffeomorphism invariance j2I] of the path integral. The 
SR generating functional is obtained by introducing the SR d'Alembertian 
□4 = ^2 + (cf. in the Appendix for d = A) and measure 

\/~—gl = X^—g2 in the classical action and by integration over the angular 
coordinates ^, (/?: 

Z\g2\ =^ (^^VXS) ■ e-^^^i'^ s{xn2+vxv)Sv^d2._ ^26) 



It is now convenient to define a new field 5" := yXS such that^ -Z^m[fi'2; S'] = 
— fj^ S'OMS'y'—g2(Px defining the complete d'Alembertian of the dilaton 

^In 1201 the original field S was preserved, leading to an additional dilaton factor X from 
the measure during spherical reduction. In the 2d action this difference can be described 
by a conformal transformation of the metric by this factor, not affecting the Hawking flux 
but other components of the EM tensor. 
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del^ 
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^ {XU2 + VXV)^ = U2 + Em, 



'X 

iyxf ux 
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2X 



(27) 
(28) 



Inserting these results into (j24j) and going back to Lorentzian spacetime dr = 
idt,WM = iWs,A — > —D,R£ — > —Rm.Es —>■ —Em the effective action of 
the dilaton model (jSj) follows^: 



96n 
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—R 
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□ 






□X 


■)■■ 


- 2 






X 



^a;. (29) 



A particular, attractive feature of the 2d dilaton model is that most of its 
effective action ()29|) can be brought into a local form by choosing a conformal 
gauge gai3 = e'^'^Va/s of the spacetime metric. The scalar curvature then 
becomes R = —2np and the endomorphism E = 00— (V0)^. It is convenient 
to represent the dilaton field in the form X = e'"^^. If one (naively) uses the 
relation 

(30) 



one obtains the local part of the effective dilaton action 

Wi[g] = ^J^{ -3co + 3ci(V0)'+pnp+6p(V0)'-6pn0 } ^/^d^x (31) 

which turns out to be identical to the one derived in ^j. However a new 
contribution 

Wni[g] = ^ { 3 N - g P<P - (V0)2] } V^d'x (32) 

appears in our present approach which cannot be brought into a local form. 
Nevertheless, we will be able to show below that the relevant expectation 



indicates Minkowski signature. 
^In the following we omit the index dil. 
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value derived from (jH^ is local after all but ill-defined. It is remarkable that 
the divergent terms do not contain the scalar curvature. 

Obviously, the form ()3H1 is not unique. Namely, by naively using relation 
(jHUjl one implicitly disregards a homogeneous solution x of the wave equation 
□x = 0. It is the aim of the next section to present a heuristic argument 
that the proper choice is indeed x = 0, and thus (jHUj) is in agreement with 
the boundary conditions which we should impose onto the scalar field S. 



3.3 Homogeneous Solution and Boundary Conditions 

The relation of a particular choice of x to the boundary conditions of the 
(massless) scalar field 5* can be seen when writing the inverse d'Alembertian 
as an integral over the Green function of 5* and applying Green's theorem 
whereby □G(x, x') = —6{x — x') and □/ = F: 



—F 

□ 



G{x,x')F{x')^^'(fx' 



dL 

oo 



- -fix) 

f{x')dr'G - Gdr'fix' 



fV^G - GVJ J y^e'-pidx'f 
dr' 



f\x')dt>G 

2M —) 



+ ... 



t'=oo 



t'=-oo 



1 ]dt' 



+ ... 



r'=2M 



(33) 



In the step to the third line the explicit form of the Schwarzschild metric 
has been used. According to (jSSl), to guarantee regularity of the boundary 
terms, G{x, x') must vanish at least linearly for r' = 2M. This is a natural 
condition for the Green functions since the manifold L is in fact a half-plane 
that is bounded by the coordinate singularity r = 2M. Further, one must 
employ some infrared regularisation to render finite the support of the Green 
functions. As the flux is measured at some finite distance r from the BH 
at some instant of time t during the quasi-static phase all boundary terms 
should vanish after that infrared regularisation. This can be realized by 
simply dropping all boundary terms'^ . But then reduces to —^F = —f 
(provided that / has at most logarithmic divergences on the horizon) and 
relation (jHUjl is fulfilled. 

It is instructive to analyze expressions like more explicitly. Unfortu- 
nately, the Green functions on a Schwarzschild spacetime cannot be obtained 



^Alternatively one could introduce a finite mass parameter or consider a manifold of 
finite size. 
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in a closed form. Nevertheless, one may construct a heuristic argument by 
considering the properties of a formal perturbation series 



G{x,x') = Go{x,x') + J Go{x,x")6D"Go{x" ,x')(fx" 
nil niii 

+ Go{x,x")6n"Go{x",x"')5n"'Go{x"',x')d^x"d^x"' + (34) 

Jl Jl 



where 

is the two-dimensional Schwarzschild d'Alembertian and 



(35) 



« = □-□„ = ^3? + * (2f a) (36) 

is the difference to Do, the flat one. 

In eqs. (ffljl . resp. (ffH|) of Appendix B the results for the second, resp. 
the third order of the perturbation series are given. 

The Green functions on the fiat half-plane which fulfill the imposed 
boundary conditions are easy to construct, singling out the appropriate eigen- 
functions. For instance, the fiat retarded Green function reads 

G^-* = e{t-t')Gl^^ (37) 

= l[e{r-r' + t' -t)-e{r-r' + t-t') 

-e{r + r' -AM + t' -t) +9{r + r' -AM + t- t')]. (38) 

In the same manner the fiat advanced and Feynman-type "causal" Green 
functions can be obtained. It should be noted that the non-local terms in 
the effective action, having the form G{x, x')F{x')^—g(P'x' with time- 
independent F = F{r), are independent of the type of Green function as can 
be verified easily. Thus only the retarded Green function (p?7j) need to be 
considered and any further ambiguity of that type disappears. 

Among the non-local expressions that appear in the effective action 
the first one has the form n~^R. The first three orders of the perturbation 
series with (jH^ and the formulas (ffUl - lfTUj) of Appendix B yield 



J G{x, x)R{x') \/ —g'd'^x 



. , 11 2M 2M2 8M3 

~— IH \ \ ... 

2 3 ^ 

-^InM -^J -ln0 = 2p-lne, (39) 



11 



e 0, which is suggested by the formal summation in the first hne. Bearing 
in mind that R = —2np = —Din (l — ^) we may argue that the per- 
turbation series produces a homogeneous solution which is just an infinite 
constant Xp = ~| Ine"- Its appearance could have been expected also as it is 
"required" to shift the absolute value of the integral in accordance with 
the boundary conditions: if one inserts r = 2M on the l.h.s. the value of the 
integral must become zero. The action of the inverse d'Alembertian on the 
half-plane on a function F = □/ thus becomes 

□-^□/(r) = /(r)-/(2M). (40) 

But this is consistent with the proper definition of a delta-function on the 
half-plane: 

6^P{x - x') := 6{r - r')6{t - t') - 5{r + r' - AM)5(t - t') (41) 

Using 5^"^ in Green's theorem (jH!^ immediately yields indeed ()4()j) . In any 
case this boundary term does not change the effective action (p?T|l because 
all homogeneous constants can be absorbed anyway by the renormalisation 
constant Ci, which, on the other hand, has no effect on the Hawking flux (c.f. 
(fS^ below). 

Finally, the first order of the second type of non-local term □"^□{/) = 
( ^^3~'' ) in the effective action should be checked: 

^"(-■)^^^^-'(w)-0-^)-- (-) 

The leading order In (2^) corresponds to the expected result as = — In r. 
To this order nothing can be said about the additional term — (l — ^) 
which should be absorbed by higher orders of the perturbation series. 



4 Hawking Flux 
4.1 Regular Part 

Having derived in Section 2 the effective action of the dilaton model af- 
ter fixing its local form (jM]) by considering boundary conditions, it is now 
straightforward to calculate the expectation values of the basic components 
by functional variation for p and X (or, more conveniently, 0). First we only 
regard the regular part oc /^[pDp + 6p(V0^) — QpOiCpl^/^dPx of (pTTjl . The 
trace © then becomes 

1 5Wreq M 
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and the pressure component (fTUI) . remembering that X = e 



5W, 



reg 



2Xy 



-9 



-1 



4M + (4M 



In 1 



2M 



(44) 

The appearance of an infinite constant In e can be traced back to the particu- 
lar boundary conditions used in this work to define the inverse d'Alembertian 
pnjl . It can be shifted to the (also infinite) regularisation constant ci, see 
next Section^°. 

As discussed in the Introduction, the basic components can be consid- 
ered as being state- independent, whereas the proper quantum state for a 
well-defined effective action has been argued to be the static (unphysical) 
Boulware state. Thus we must rely on the conservation equation (jHl) to 
compute the remaining components of the EM tensor in the Unruh state. 
A necessary condition for regularity of the EM tensor at the horizon was 
Qu = 0. The other constant Ku is determined by the condition that the 
incoming fiux vanishes 

fl_ 2M) 

{U IT++I U) = {{T\) - {T-'J - 2 {T\)) 



{{T)-2{T\)-2{T^*t)) 



r— >oo -''■[7 
^ M2 



Krr irMmiT) + {l-m)d^X{T\) 



dr' 



Qu 

2M2 



0. 



(45) 



Inserting (j4Sj) . ()44|) and fixing the dilaton field to its standard form X = r 
in spherical coordinates we obtain 



u 



2M 



M 



2(r 



i\2 



(T), + (r' - 2M) (T%>^ 



dr' 



1 



7687r 



(46) 



By the CP equations (0),® the total fiux F^^g from the regular part of the 
effective action through a large spherical shell surrounding the BH is then 
given by 

^ -{u\r\\u),--^- — - 



which precisely coincides with the result obtained in refs 
related fiux in 4d through a sphere of size 47rr^ becomes 

1 



Uni EE]. 



(UlT'-tlU), 



30727r2M2r2' 



(47) 
The 

(48) 



^"For the moment we keep this notation until we consider the contribution of the cor- 
responding divergent term in the effective action. 
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As a consequence the BH behaves as a black body at Hawking tempera- 
ture with the radiation flux according to Stefan-Boltzmann's law [T^ I16j. 
It is appropriate, though, to emphasize at this point the drawback of this 
solution to the energy flow problem at the horizon the energy density 
and other components of the EM tensor exhibit a logarithmic singularity in 
global coordinates! As argued in [12] this singularity (being integrable) is 
not in contradiction with the finiteness of the total flux. Actually, we also 

obtain the outgoing flux in light-cone coordinates T with that singularity 

(eq. (97) of [ini) in our present approach 



(l - 

= 768vrMV ^ + + ^^^^ + 



2M , , 
In 1 ] - \ne 



(49) 



4.2 Divergent Terms 

The first divergent term oc Cq = e"^ in (|^ is a pure UV divergence and has 
the form of a cosmological constant: 

Wco = ^ [ coV^d'x = ^ / e'^'d'x (50) 

It contributes to the trace of the EM tensor (T) = ^ and hence also to the 

^ I Co An 

asymptotic flux {e has dimension M^) 

^ " L - 16vr " ^^^^ 

This divergent contribution to the EM tensor can be interpreted as an infinite 
vacuum energy because it even appears in the case of flat spacetime M = 
R = 0. A renormalized EM tensor can thus be defined by subtracting the 
flat spacetime value (with flat metric r]) (T^p)^^^ := {Taf^)g^ - {Tais)^ or by 
simply setting cq = 0. 

Further, we had an IR-UV divergence cx Ci = ln(£:/T) in (jSH), contribut- 
ing only to the pressure component 

Although the asymptotic behavior of the EM tensor is unaffected by ci be- 
cause 

poo 

K^^ = MM (r - 2M) (T%)^ dr = (53) 
J2M '^^ 
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it produces infinite contributions to the EM tensor in higher orders in r. 
Comparing with (jl^ we observe that Ci appears in the same place as the 
homogeneous solution Xp = — | Ine needed to fulfill the boundary conditions 
(cf. Section 3.3) by (jlO)). After all, any constant homogeneous solution Xp 
leaves the asymptotic flux invariant and can be shifted to the regularisation 
constant Ci. 

The situation clearly is different if the homogeneous solution is a function. 
If, for instance, one chooses as in ref. HI 



Xp 



2M 



1 + ln 



2M 



- 1 



Ine 



(54) 



in order to eliminate the logarithmic singularity of the flux at the horizon, 
this would mean that the corresponding non-local term in the effective action 
had the form (cf. 



G{x, x')R{r')(fx' 



2M 



-1+ln 



2M 



1 - 



2M 



-2 Ine, (55) 



and the asymptotic Hawking flux not only is affected, but even would become 
negative because Ku > 0: 



Ku = Ku + M' 



(r'-2M)(r' 



2M 



^^Xpdr' 



Ku + 



1 



1287r 7687r 



(56) 



This is the result obtained in ref. eq. (16). Also the other components 
of the EM tensor calculated with xp in ()54|1 can be verifled to agree with 
those of that work. As observed by the authors of themselves, then the 
weak energy condition is clearly violated in the asymptotic region. We do 
not believe that this serious consequence of the choice justifles its aim 
to eliminate the logarithmic singularity at the horizon [10^. 

Finally, we have a non-local divergent part ()32|) of the effective action oc C2 
whose contribution to the basic components can be localized in a conformal 
gauge 



(T) 



1 6W, 



nl 



1 



C2 



^ ln(-Te-^^Do) , 
E = Ey'—gd X 



□ 



^^^^ ^ np[p - p{2M)] ^ M[ln(l-g^)-lne] 
47r □ An Anr^ 



1 



5W, 



nl 



1 



C2 



2X. 



C2 ^ 
-Ey 

66 □ 



-gcPx 



(57) 
(58) 



[□ + 200 + 2V0V] C2 ^ [□ + 200 + 2V0V]c2[p - p(2M)] 

-E = 



87rr2 



□ 



87rr2 
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where we have used the relation (after having varied the effective action!) 



, AM-r 2M-r 2M R 
E = 00- (V0)^ = -y^ ^ = _ = __ = Dp. (59) 

()58|) contains an ill-defined expression In(Tn) that cannot be treated further. 
Even if C2 were only a constant, the pressure component would change by 
a term {T%)^^ = g^{4M + (4M - r)[ln (l - ^) - \ne]} identical to the 
original expression (jH}. In that case the asymptotic flux would be affected 
too: Freg — > Freg " (1 — 3c2). Bccausc Freg is supposcd to be the correct result 
this constant had to be zero. 

The appearance of IR divergences within the covariant perturbation the- 
ory (in contrast to the local Seeley-DeWitt expansion used in former deriva- 
tions fSJ ) is not surprising as it allows for infinitely large values of the 
eigentime r, corresponding to zero modes. As long as a more detailed anal- 
ysis of these terms does not exist we can only highlight the existence of such 
terms while assuming that their (infinite) contribution to the EM tensor 
might be canceled by some mechanism that has not been considered yet and 
may necessitate the inclusion of (even arbitrary ?) higher orders in covariant 
perturbation theory. 

5 Conclusions and Outlook 

The Hawking flux from a spherical Black Hole has been reconsidered, whereby 
we followed the line of solving the EM conservation law, as proposed by 
Christensen and Fulling many years ago ^21- We have linked the formalism of 
d = Ato the "non-conservation equation" in the effective d = 2 dilaton theory 
where as in 4d gravity, beside two constants, only two "basic components" 
are a necessary input, the 2d trace anomaly (T)2 of the EM tensor and the 
"dilaton anomaly" {T^e)2 to be interpreted also as (part of) the pressure 
component (T^e^^ = {T'f^)^ in = 4. 

Whereas the computation of (T)2 is known for a long time to fit perfectly 
into the heat kernel formalism, in its single previous determination of refs. 
[ISl EE] {^^9)2 '^^^y ^ ^^^^ force argument had been made accessible 
to that technique. In our present work we replace both derivations by a new 
application of the "covariant perturbation theory" , introduced by Barvinsky 
and Vilkovisky fr7\ ITH] , which allows the direct determination of the effec- 
tive action also in the presence of dilaton fields. Again the heat kernel is 
used, albeit in a slightly different manner. The non-local form of this ac- 
tion implies an important dependence on boundary conditions for the Green 
functions of the scalar field in the given background. We argue that the 
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effective action only has a consistent interpretation when for that quantity 
the (otherwise unphysical) Boulware state determines the boundary condi- 
tion at infinity (infrared regularisation) . This is supported by a study of the 
Green functions where - in a formal perturbation series starting from flat 
spacetime - a natural boundary condition at the horizon directly leads to a 
constant homogeneous solution. The latter can be absorbed in a renormali- 
sation constant which does not affect the asymptotic flux. Generically, the 
choice of different homogeneous solutions represents the ambiguity inherent 
in the non-local effective action. 

From the regular part of that effective action in the Boulware state the 
Hawking flux is not derived directly, but only the non-radiative basic com- 
ponents (T)2 and (T^q)^ which are, indeed, independent of the assumed 
(Boulware, Hartle-Hawking, Unruh) quantum state [20^ . Inserting them into 
the EM (non-) conservation equation the correct Hawking flux in the Unruh 
state (in agreement with the black body derivation) is obtained, without pro- 
ceeding through one mathematically questionable step which in refs. |13| IT^ 
produced the same physically reasonable result. 

We also clarify the consequences of the ambiguity brought about by dif- 
ferent other choices of that homogeneous solution which are not constants. 
They may even yield a negative flux ^O]- A choice like the one in that 
work, as noted by the authors themselves, violates the weak energy condi- 
tion which we consider a more serious defect than the (logarithmic, hence 
integrable) divergence of some of the components of the EM tensor at the 
horizon QEj • 

Beside the above-mentioned attractive features of the covariant perturba- 
tion theory, it reveals the existence of three new divergent terms (up to this 
order in curvature) which require further investigation. The first one oc Cq is 
interpreted as an infinite vacuum energy and its contribution to the flux can 
be removed by common renormalisation arguments. The second one oc Ci 
seems to be related to the boundary conditions as it contributes another log- 
arithmic (UV) divergence at the horizon, as compared to the one produced by 
the homogeneous solution Xpi leaving a logarithmic IR divergence. Although 
we could show the independence flS^ of the asymptotic flux of Ci as well as of 
any constant Xpi the presence of an infinite constant in higher order terms in 
r could not be avoided (at least, when renormalizing this constant to zero one 
is left with the logarithmic divergence at the horizon as in jEllISl)- Finally, 
our approach yields a completely new term for the effective action that 
cannot be localized even in conformal gauge and contains an IR divergence 
C2 which is coupled to an expression In □. Although it could be shown that 
its contribution to the EM tensor can be given in a local form, the latter 
could not be evaluated further due to its ill-definedness. To sum it up, the 
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covariant perturbation theory (by its non-local character) produced new IR 
divergent terms that were absent in the local Seeley-DeWitt expansion used 
in former derivations ^HlEj- A proper renormalisation of these terms would 
require a detailed analysis of this method (e.g. considering higher orders in 
curvature or partial summations) which might be the issue of future work. 

Beyond the application to the present problem the extension of covariant 
perturbation theory may set an example for the use of that technique in 
theories where matter is also coupled non-minimally to a scalar field (such 
as scalar-tensor theories). Also applications to a non-perturbative approach, 
where the geometry is integrated out, seem promising [221 ES]- 

The improved understanding of the situation for a static spherically sym- 
metric Black Hole also seems a good basis for extending our result to include 
e.g. "grey factors" (ref. 50 in |F) by admitting a non- vanishing classical 
scalar background field, back reaction upon the metric, and also the consid- 
eration of higher orders in the covariant perturbation series. 
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A Spherical Reduction 



The spherical reduction procedure is the basis of the two-dimensional dila- 
ton model considered in this work. We reproduce it shortly, because some 
intermediate formulas are important for the main text. Here we consider a 
more general reduction from a ci- dimensional spacetime M with spherically- 
symmetric metric 



ds^ = gafjdx^dx^ — ^'^{x'^)g^\dx'^dx^ 



(60) 



to a two-dimensional Lorentz submanifold L, spanned by the coordinates 
(e.g. t, r), by reducing out a. {d — 2)-sphere S''^"^ with coordinates x'^ = 
6, If ... . The dilaton field is defined as X = <l>^, $ is more convenient for 
calculations. gaf3 is the induced metric on the L and g^x the one on S'^^'^. 
We work in a vielbein basis in which the line-element can be written as 



ds^ = riabe" ® 



(61) 



The form a vielbein basis on M. One can define a vielbein basis on L and 
S'^~'^ which we denote by and e*, respectively^^. They are related to the 

e'^ by 



e = 
Ea , Ei 



= ^-'Ei. 



(62) 
(63) 



Further, a (Levi-Civita) spin-connection uj'^s on M induces connections on 
the submanifolds: 

Co^b = uj^b , ^'j = uj'j- (64) 
The connection on M, determined by vanishing torsion on M and L, is given 

by 

With ()65|) the Riemann tensor on M can be expressed by geometrical objects 
on L and S'^~'^. For the scalar curvature one has 



(65) 



R 



M 



+ {d - 2){d - d.) [Eb^] iP^ 



R' 



EbE% 
{d -2){d- 3) 



E'<^]Co\ E 



$2 



1 + Vfe^V^ 



d-2 



□l^. (66) 



^^We mark geometric objects belonging to L ot S'^ ^ by a tilde; tensorial objects on the 
submanifolds are distinguished easily by the different indices used. 
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In the last line we have inserted the (constant) scalar curvature of the {d — 2)- 
sphere = {d — 2) {d — 3) (the Riemann tensor on the (rf— 2)-sphere is given 
by R^j = e* A In this form all quantities on the r.h.s. live on L as it 
should be. By reduction from d = 4 the Ricci tensor and the scalar curvature 
become 

RM ^ jj.„2i±(|£)!_4™=fi^-l+(Z^-2™. (68) 

In the last equality we have returned to the dilaton field X = $2 Note 
that the ^-term in the scalar curvature becomes a surface term in the 
action when it is multiplied by the SR measure y/—gM = XyJ—gi. For 
a d'Alembertian on the c?- dimensional manifold M acting on a scalar field 
S{x°') which depends only on the coordinates of L one gets 

US = r]'''VrEsS = t^VaEbS + r]'^ViEjS 

= DS- r,'^uj'^,{E,)EaS = nS~ V^'^KS 

^□5+(.-2)X!^.n^ + l^™^. (69) 
$ 2 X 

If M is the four-dimensional Schwarzschild spacetime and the gauge of the 
dilaton is fixed as X = this simplifies to 

DS = DS-'^(^1- drS. (70) 

B Green Function Perturbation - Higher Or- 
ders 

In this Appendix the second and third order of the Green function pertur- 
bation series are adapted to the case of a 2d Schwarzschild spacetime 
with perturbing d'Alembertian (jSE)), taking into consideration the particu- 
lar boundary conditions on the half-plane and the time-independence of the 
involved integrals. The second order of the perturbation series then reads 
{d" = dr") 

r" 

[d"Go{x, x")]g{r")d"Goix", x')d^x", (71) 
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where one r"-derivative has been partially integrated and 

. s 2M 

g{r) := —. (72) 

On a two-dimensional Schwarzschild spacetime d^g(r) = —R{r). A useful 
identity can be derived, introducing flat light-cone derivatives d± = dt ± dr , 
□o = d+d- = d-d+. Because of the time-independence one can effectively 
set Do = -d^. and d+ = -d- = dr, {g" = g{r"),Go{x,x') = G^cx'^^oGxx' = 

. ^ - / + 2.",a"G..,,a"G..} ,73) 

By ()73|1 the second order of the perturbation series ()H4|1 can be written in 
the compact form 



^i^[g(r)+gir')]Go{x,x') + J Go(a;, x")Go(x", x')i?(r")rfV| . (74) 
In a similar manner the third order is computed, 

/// /•/// 

/// /•/// / <•// 



(75) 
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where again the identity (|7S|) has been used. Finally, we write down the third 
order of the perturbation series in a compact form, replacing drg'^{r) = 2R{r): 

I I [g'{r) + gir)gir') +2g'ir')]Goix,x') 



+ J Go{x, x")R{r") [2 ■ dr" + g{r')]Go{x" , x')<fx" 

/III 
Go{x"',x')R{r"') 

[g{r) + g{r"')]Go{x^ x'") + / Go{x^ x")Go{x'\ x"')R{r")(fx" ) cfx'" 



(76) 



C Non- Conservation Equation 



Diffeomorphism invariance of the generating functional implies a non-conser- 
vation equation for the renormalized expectation value of the EM tensor when 
a dilaton field is present, resulting from the non-minimal coupling of the 
scalar field in c? = 2 (cf. Q). A 2d diffeomorphism transformation S^g"'^ = 
-(V"^^ + V^r), = r^a^, = ^"d^S applied to the effective action 
W[g] = —ilnZ[g] yields the quantum non-conservation equation which has 
not been checked in the previous literature: 



6i:W[g] 
-iM 



i5i:Z[g\ 



iM 



m 

-iM 



m Z[g] 



Z[9\ J JL 



-9r 



e d y 



'-aid 



+Cd^x\^ [{VSf - m'S'] d'y 



iM 



Z[g] 
M 



V{^gS) e^'^ l'±6^Sd'y 



I V{^gS) e^- yl^'r {v^T„, + ^ [(V^)^ - m'S^] 



lim da6{x - y) \y/^d y 



-g 



= lim da6{x - y) >V^d y 

-g x^y 



(77) 
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In the step from the second to the fifth fine a partial integration in the path 
integral has been performed. The divergence 5(0) represents the infinite 
zero-point energy of the quantized scalar field S. Defining the renormalized 
EM tensor {Tap)^^^ := (Tq^j) + ^= Ymix^y 5{x — y) we obtain indeed (jHI), as 
proposed. 



D Heat Kernel Integrals 



The leading divergent (T — > oo,£: ^ 0) terms of the expansions ()19II23|) . ap- 
pearing in the trace of the heat kernel (fTTj) can be computed analytically. For 
instance ()2()j) only has an UV divergence, hence we cutoff the r-integration 
at the lower boundary by e > 0. As a first step we expand the regularized 
integrand in powers of s, carry out the differentiation for s and then set it 
zero: 



d 



1 



ds 1 r(s) 
d 

ds 



dr 



s=Q 



|[s + s^^E + Ois^)] j^[r-^ + Inr + 0(s2)]rfr| 



s=0 
1 



T-'dT = -. {r 



The most problematic term is ()21|). Introduction of a cutoff at infinity and 
differentiation for s leads to (using (|T8|) : 



d 



ds Vis) 



T'f{-TA)dT 



s=0 JO JO 



1 g-ra(l-a){-A) _ 

a(l-a)(-A) 



da 



-ra(l-a)(-A)^^^^ 



1 e(-^-i)^ - 1 



-A 



dz (79) 



This integral cannot be solved analytically, however, it can be compared to 



^0 



' 4 drdz 



-A 



1 e(-i)^ - 



z-\ 



In 



-dz 



-A 



Tf-A) 



+ lE 



(80) 



This is twice the result conjectured already at the r.h.s. of (PT|) apart from 



an additive constant Hence, it remains to show that the difference 

between //2 and of the original expression (l.h.s. of (I2H)) converges to (one 
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half of) that constant for large regularisation parameter T. First we perform 
a substitution of the integration variable z ^ in (IHUj): 



2 / / 



^ r r z-e^^'-^y-^dzdT (81) 



This difference then becomes 

4 (^ _ l)dTdz = —- / -dz 



e ^ 

^0 -^Jo z + 1 



A Jo ^ + 1 A A 

which proves (j2T| to be the correct result for large T. We note that the limit 
limr^oo e^^^~^^ 2 = could be performed because < 1 for all z except 
for z = 1. But at that value the integrand vanishes altogether. 

The computation of the remaining expressions is tedious but straightfor- 
ward. 
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